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Abstract 

The ATP-driven Plasma Membrane Calcium pump or Ca 2+ -ATPase (PMCA) is char- 
acterized by a high affinity to calcium and a low transport rate compared to other 
transmembrane calcium transport proteins. It plays a crucial role for calcium extrusion 
from cells. Calmodulin is an intracellular calcium buffering protein which is capable 
in its Ca 2+ liganded form of stimulating the PMCA by increasing both the affinity to 
calcium and the maximum calcium transport rate. We introduce a new model of this 
stimulation process and derive analytical expressions for experimental observables in 
order to determine the model parameters on the basis of specific experiments. We fur- 
thermore develop a model for the pumping activity. The pumping description resolves 
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the seeming contradiction of the Ca 2+ :ATP stoichiometry of 1:1 during a translocation 
step and the observation that the pump binds two calcium ions at the intracellular 
site. The combination of the calcium pumping and the stimulation model correctly 
describes PMCA function. We find that the processes of calmodulin-calcium complex 
attachment to the pump and of stimulation have to be separated. Other PMCA prop- 
erties are discussed in the framework of the model. The presented model can serve 
as a tool for calcium dynamics simulations and provides the possibility to characterize 
different pump isoforms by different type-specific parameter sets. 

Keywords: Plasma Membrane Calcium pump, Plasma Membrane Ca 2+ -ATPase, calmod- 
ulin, stimulation, relaxation, pumping activity, theoretical model, parameter 



Introduction 
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asma Membrane Calcium pump with a high calcium affinity {K\n < 0.5 — 1 fxM 
2|) and a low transport rate (~ 30 Hz, J] and private communication) is an important 
component for the maintenance of calcium homeostasis in cells. By using the energy stored 
in ATP the PMCA transports intracellular calcium ions out of the cell. It has been found 
in all mammalian eels 5j, where the expression level does not exceed 0.1% of the total 
membrane protein . An exception is the brain where this value is up to 10 times 

higher than in non-excitable cells jf|. 

The four different pump isoforms are encoded by four independent genes, which are 
indicated by numbers 1-4. The diversity of pump forms is further increased by alternative 
mRNA splicing variants, characterized by small letters. There exist more than 26 transcripts 
which differ in their regulatory properties, for instance in their affinity to calmodulin, and 
which are distributed in a tissue specific manner Q, [J. Referring to the available data we 
will investigate hPMCA2b and hPMCA4b, where h stands for human j^|. 

Calmodulin is an intracellular calcium sensor protein with four relatively high affinity 
Ca 2+ binding sites (K a = 1 — 10 /iM" 1 at low ionic strength It belongs to the mobile 

proteins of the EF-hand family with a helix-loop-helix conformation Q|. Two calcium ions 
are bound at the N-terminal as well as at the highly homologous C-terminal domain, each 
of them formed by two EF-hands M. The domains are connected by a flexible linker - an 

n 

a-helix Because of its calcium binding capacity calmodulin becomes relevant for the 
spatial propagation of calcium signals within the cytoplasm. Additionally the fully liganded 
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calmodulin-calcium complex is responsible for Ca 2+ -dependent regulation of the activities 
of a vast array of different target proteins, including enzymes, ion pumps and channels 
Q|. Amon g th ose with high affinity (i-Q = 5 nM) to the calmodulin-calcium complex is 
the PMC A |l.T |. The free calmodulin-binding domain of the PMC A also interacts with 
the ATP-binding site of the pump and acts as an inhibitor of ATP-driven pump function 
Q The detachment of that autoinhibitory domain segment after binding of the complex 
causes a stimulation of the pump function bv increasing both the affinity to calcium and 
the maximum turnover rate of calcium P, 1^,0, Q|. The stimulation and the relaxation to 
the initial unstimulated state happens on a time scale of minutes and enables the pump to 
display a memory of previous calcium transients j^|. 

Caride et al. have published a stimulation model making use of measurements of the 

n 

dynamical stimulation behavior [2J. We extend the reaction scheme and aim to improve the 
model results in two ways: Firstly, we will deduce the stimulation parameters from mea- 
sured data. Secondly, we will include a saturation of the stimulation rate for high calcium 
concentrations. Indeed, measured stimulation constants call for a limitation by a maximum 
value at high calcium concentrations. After completion of this work we were made aware 
of the recently published article in which Penheiter et al. propose a new stimulation model 
in conjunction with fluorescence measurements Q|. Although Penheiter 's model includes 
saturation, the relaxation of the PMCA to resting state is not considered. We separate the 
stimulation as well as the relaxation into two steps and introduce rate limiting reactions in 
our model. An analytical approach enables us to deduce the required rate constants from 
measurements. The inferred system of coupled nonlinear differential equations is solved 
numerically. The presented new model is in agreement with the investigated experimental 
data of the stimulation process by calmodulin. Furthermore we predict the stimulation 
behavior beyond the available data. 

Besides the stimulation model we also consider the calcium pumping activity. This 
joined model of stimulation and pumping may serve as a tool for simulations in a wide 
variety of systems such as single cells or tissues. For that purpose we assume that all 
isoforms are functionally similar and characterize the different isoforms by type specific 
sets of parameters. In this way we differentiate between universal system independent 
parameters which are characteristic for each pump type and specific parameters such as 
the PMCA protein expression level, which has to be adapted to each experimental setup. 



In relation to the available experiments we can specify parameter sets for the h2b and h4b 
isoform which turn out to be in a biologically reasonable range. 

We begin by introducing the stimulation model reactions. We determine the model 
parameters and calculate stimulation dynamics. In the second part of this paper we focus 
on the derivation of an expression for the calcium pumping activity. Finally the stimulation 
and the pumping model together serve to fully describe PMCA function. 



Model 

Stimulation model 

We base the stimulation model on a system of rate reactions, which describe the calcium 
and calmodulin dependent transition from unstimulated to stimulated pump form and vice 
versa. The stimulation dynamics of an ensemble of PMCA's is characterized by the follow- 
ing reactions: 



P x ^ P* + X 4 , (4) 

fe_4 

P* ^ P. (5) 



We use the following abbreviations: P denotes free PMCA, CaM-Ca4 = X4, P-CaM-Ca4 = 
Px, P*-CaM-Ca4 = P x , and P* free stimulated pump. P and Px are the unstimulated 
pumps and the asterisk denotes stimulated pumps. M and Ca 2+ are free intracellular 
calmodulin and calcium, respectively. Note that in this article Roman style symbols refer 
to elements or proteins whereas italic style symbols denote concentrations or fractions of 
the respective element or protein. The fc's are the rate constants. With the Law of Mass 
Action the model CQ)-© is rewritten as a set of differential equations for concentrations 
(see equations (jf"6|) - (fl"8)) in the Appendix for the stimulation model). 



M + 4Ca 2+ 


feTi 


X4, 


(1) 


P + X 4 




Px, 


(2) 




fcT 2 




Px 


hi 


Px, 


(3) 



Reaction Q describes the binding of calcium to ca^ 



modulin. The cooperative binding 
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is simplified in this reaction by 



of Ca 2+ within each domain of the calmodulin protein 
the assumption of highly cooperative binding of all four calcium ions. Based on the realistic 
assumption that the binding of calcium occurs faster than the PMCA stimulation, which 
happens on a time scale of minutes, we treat this binding to be in quasi-steady-state. 
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The assumption of irreversible stimulation Jjj} and relaxation © implies that the pump 
can only be stimulated when X4 has attached to the pump, whereas recovery back to the 
unstimulated pump form only occurs when X4 has detached. It will be seen later on that 
the irreversible step J2J and the detachment step (@J limit the stimulation and relaxation 
rate, i.e. both cannot exceed the rates k^ and k^, respectively. Note that (@J does not 
necessarily mean that the complex detaches as a whole. In principle the calcium ions could 
also detach from first without any alterations of the model because of the quasi-steady- 
state approximation in (£Q). Thus, the essential statement of (@J is that calmodulin (with 
or without calcium) has to detach in order to destimulate the pump. 

Note that we introduced different pathways for the stimulation (reactions (J2J) and (PJ) 
and the relaxation process (reactions (HJ and (JHl). This enables us to describe the temporal 
behavior of these two processes independently of each other, which is rather important in 
order to reflect isoform specific stimulation and relaxation dynamics. 



Stimulation model parameter determination We aim to relate the rate constants 
ki (i = ±1, ±2, 3, ±4, 5) and the stimulation and relaxation constants k^^ and k^ ax 
measured in experiment by Caride et al. jj. Note that in [^| these constants are denoted 
by k ac t and k- mac t, respectively. Starting from the stimulation model we deduce analytical 
expressions for the exponential growth constants &; s tim and /c rc iax- 

Stimulation Suppose only unstimulated pump form P to be present with the subse- 
quent addition of calmodulin. In the very beginning the stimulation model can be reduced 
to the rate equations |[2j) and (0- The binding of calcium to calmodulin (reaction 
£Q)) is assumed to be in quasi-steady-state, therefore X4 = CaM-Ca^ = K ■ M ■ Ca 4 , 
with K = k\jk^\. During the beginning of the stimulation the finite calmodulin concen- 
tration imposes no restrictions on the dynamics, since only a minor fraction of free X4 
is bound. Therefore, the conservation of the calmodulin concentration is neglected and 
M is considered to be constant. With this assumption and the pump mass conservation 
Po = P + Px + Px (note that P* is not involved here) reactions @ and (JBJ can be rewrit- 
ten as a system of three linear homogeneous differential equations of first order (the full 
nonlinear system is shown in equations pfi [l -l|18 )) in the Appendix). The system is solved 
via the ansatz P = A\ ■ exp(—k st i m t), Px = A2 ■ exp(—k s n m t) and P x = ^3 • exp(—k s f im t). 
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The characteristic equation leads to three solutions for k st - mi : k} ti = and 



, (2/3) fc-2 + h + k 2 X 4 fk^ 2 + h + k 2 X 4 \ 2 

fc stim = 2 V I 2 J ~ /C2/C3X4 - (6) 

With the boundary conditions P(0) = P , P(oo) = 0, P x (0) = 0, Px(oo) = 0, P£(0) = 0, 
PJf (oo) = Po and the pump mass conservation Po = P(i) + Px(t) + PJK*) ^ the exact 
solution for the stimulated pump form becomes 



P*x(t) 



V^stim / \ ^stim/ 



(7) 



During the first phase of stimulation the pump exhibits a single exponential behavior 
as discussed in Q|. The exact solution for P x (t) in equation Q can be simplified to a 

(2) (3) 

single exponential expression in two cases. With A; stim >> fcg tim equation Q becomes 

Pl(i)=P (l-e-^l*). (8) 

(2) (3) 

Formally, also & stim << fegtim yields a single exponential behavior of the same form but 
with fcg tim in the exponent. However, k^ t ( m > & stim holds true (see equation ©). 

(2) (3) 

Assuming /Cg tim >> ^ stim the formation of the stimulated pump form during the first 

(3) 

phase of stimulation depends on fc stim only. This approximation is supported by the fact 

that fcgy m converges at high calcium concentrations in the same manner as k^^ does (see 

figure Pi . In contrast, k^ t ( m diverges for high calcium concentrations (see Appendix). & stim figure ^ 

can be interpreted as corresponding to the single exponential fit constant k^^. This enables 

us to determine the stimulation parameters (shown in figure Q). 

In an analogous fashion to the stimulation case we deduce an analytical expression like 
equation (JHJ) for the relaxation scenario. The exponent of this single exponential expression 
can be related to the experimentally measured relaxation constant (see figure El calculation 
not shown). figure [21 

A fit of the derived expressions for the introduced stimulation parameters to the ex- 
perimental data restricts their values. As a first approach we use Hill equations (see Ap- 
pendix for the stimulation case). The maximum stimulation and relaxation fit constant 
can be identified with k% and k±, respectively, which limits the stimulation and relaxation 
rate. In contrast, the half maximum concentration of the Hill fit is not sufficient to deter- 
mine the remaining free parameters. Even with the knowledge of the dissociation constant 
K2, K2 = k_2/k2 = 0.5 nM for the h2b isoform and K2 = 5 nM for the h4b isoform 



t Il5l. 1 1 Si. lift l2nl|. the system remains under-determined because K is not precisely 
known. A lower boundary for K comes from the requirement of positive reaction rates &2 
and /c_2 and from the value of K2. An upper boundary for K is imposed by the choice 
of the value for K4 = k i /k_ 4: . The parameters resulting from the Hill equation fits are 
reported in table table Q 

Attention should be drawn to the fact that Ca 2+ is bound by polar amino acids of 
calmodulin and therefore the affinity represented by K strongly depends on the ionic 
strength of the experimental solution We assume the same value of K for all ex- 

periments, i.e. the same conditions in all experiments. Since the value of K is not exactly 
known for the experiments under consideration, we assume three different values and quote 
the respective parameter sets in order to show how other parameters depend on K. The 
values used for K are 0.1, 1 and 10 /uM -4 which, according to Linse et al, correspond to 
an experimental solution with a KC1 concentration of 18, 20 and 26 mM. 

Using © and the parameter values determined by the Hill equation fit we are able to 
calculate the calmodulin dependence of fc st ; m without further assumptions. The result in 
fie nre □ (dotted line) is in qualitative agreement with the experimental data of figure 3B 
in Q] (A: act ). The measurements have been carried out at a calcium concentration of 10 
/uM with 0.005-0.065 //M calmodulin. The experimental data could be fitted by a linear 
equation. The saturation suggested in figure ^ occurs at higher calmodulin concentrations 
only. A linear fit is also suggested in figure 4A in for a range of 0.005-1.25 fiM calmodulin 
at a constant calcium concentration of 1 fiM. However, the linearity of the fit relies on a 
single data point and we consider this result to be in agreement with the saturation predicted 
by our model. We quote k st i m at a calcium concentration of 0.8 //M for the h2b isoform and 
at 1 ,uM for h4b. At both cases the stimulation constant saturates already at low calcium 
concentrations, i.e. at 0.2 fiM for h2b and at « 0.4 fiM for h4b. The results from the Hill 
equation fit also allow us to predict the calmodulin dependence of the relaxation constant 
fcrelax, which has not been measured so far (see figure EJ. 

We have seen that it is not possible to reliably determine all stimulation parameters by 
a fit on the basis of Hill equations. We therefore use a Metropolis algorithm, which is based 
on a least-squares fit routine, to incorporate three sets of available data for the stimulation 
constant, the relaxation constant and, in addition, the steady-state pumping activity, which 
will be discussed subsequently. The medians and the confidence intervals in table [21 and 0] 
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are calculated from a Bootstrap method, i.e. random data sets of the same size are drawn 
from the original data set from which, using the Metropolis algorithm, the fit parameters 
are calculated. Repeating this procedure 10.000 times enables us to calculate medians and 
confidence intervals for the fit parameters. The medians of the parameters in table 01 are 
determined similarly but for different values of K . 

Note that we have derived equation © only by considering the beginning of the stimu- 
lation. An alternative approach to determine the parameters would be the fit of the exact 
model equations to the time course of inorganic phosphate Phj release (whereas sub-i refers 
to inorganic). Phj emerges by hydrolyzing ATP during pumping (see section "Results" and 

a a. 

However, this would not improve the accuracy of the parameter determination 
since these data contain information not only about the stimulation but also about the 
pumping of the PMCA, which cannot be disentangled. Only a measurement of the exact 
stimulation and relaxation behavior would provide relevant additional information in order 
to get a more reliable choice in table EJ 

Stimulation dynamics Based on the stimulation model reactions (EQl-(© we can simu- 
late the dynamics of the unstimulated and stimulated PMCA states. With the obtained 
isoform specific parameters differences between the isoforms can be discussed. The fact 
that the value of for the h2b isoform is more than twice that of the h4b isoform reflects 
faster stimulation, whereas relaxation is faster for the h4b isoform since k% Ab > k% 2b . The 
stimulation and relaxation dynamics of both isoforms are shown in figure 01 In graph 01 figure 01 
(a) we start with no stimulated pump (i.e. / s tim(^o = 0) = 0) and expose the pump to a 
high calcium concentration. In the relaxation graph 01 (b) all pumps are stimulated (i.e. 
/stim(io = 0) = 1) and the calcium concentration is decreased to 0.1 /iM. 

Pumping model 

Up to this point we have described the temporal transition between unstimulated and stim- 
ulated pump forms. In this section we will mathematically characterize how the PMCA 
transports calcium ions across the membrane depending on the intracellular calcium concen- 
tration. We will consider the pump cycle to be in quasi-steady-state, i.e. the pump activity 
reacts immediately to the intracellular calcium concentration. This simplification is justi- 
fied since the pump activity adapts to different calcium concentrations within milliseconds 
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table El 
table 01 



while the stimulation and relaxation happens on a time scale of minutes. 



Pumping activity In contrast to the Ca 2+ pump of the sarcoplasmic reticulum (SERCA) 
the stoichiometry between transported Ca 2+ and hydrolyzed ATP of the PMCA is most 
likely 1:1 IJEIEIQ. The SERCA pump works with a calcium:ATP ratio of 2:1 QQ. 
Hence, the transport step of one calcium ion through the plasma membrane becomes 

Pin 7 Pout + Ca 2+ out . (9) 

Pin stands for a state of the pump in which one or two calcium ions are bound and which 
can perform the calcium transport process by a transformational change. P out refers to the 
pump state after the translocation step. We do not consider the transformational change 
from P out to P; n , which closes the pump cycle. Since ATP is assumed to be sufficiently 
available in cells we treat the ATP concentration as thermodynamical bath. The calcium 
transport rate J in the outward direction is given by 

J = n • Pin- (10) 

Carafoli pointed out that the activation of the ATPase by calcium and the saturation 
can best be described by a Hill equation with a Hill coefficient of two P|. This may be 
related to the binding of 2 calcium ions at the intracellular binding site of the pump even if 
only one calcium ion is translocated. In the following this is incorporated into the model: 

P + Ca 2+ 3P-Ca, (11) 
P-Ca+Ca 2+ 3 P-Ca 2 . (12) 

r-3 

The pump state Pi n is assumed to comprise the pumps with one and two calcium ions, 

i.e. Pm = P-Ca + P-Ca2. This implies that the pump with one calcium ion, P-Ca, is able 

to perform the transport process which is not necessarily the case but the most general 

assumption (see also comments after equation ifTljl ). 

In the quasi-steady-state approximation of equations (|TT|) and lfT2*)l we find an expression 

for Pin- Using this expression the calcium transport rate (equation <fTfl|0 becomes 

= R 3 Ca + Ca 2 

1 R 2 Rs + R 3 Ca + Ca 2 ' [ ' 

with i?2 = f-2/f*2 and R3 = r^/r^ and the pump mass conservation 1 = P+P-Ca+P-Ca2- 

Note that we have divided the pump mass conservation by Po so that the different pump 
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states become fractions instead of concentrations. The sum of the fractions of all the three 
pump states is of course 1. 

Fitting equation (|lHp to the experimentally measured transport rate in the absence of 
calmodulin yields R2 = 38 fiM and -R3 = 0.011 /jM. for the h2b isoform (see figure^}. These 
values suggest a cooperative binding of both calcium ions, i.e. the binding of the first ion 
is slow but of the second fast. In the model this means T2 — > and — > 00 while keeping 
f2 ■ 7~3 constant, in which case R2 — > 00 and R3 — > while R2R3 is constant. Experiments 

same sigmoidal behavior 

of the transport rate as found by equation lfT3|) in figure QJ supporting the assumption of 
cooperative binding. Thus (fT3)) is simplified to 

1- Cq2 - r Ca " 
- ri r 2 r 3 + Ca 2 ~ Jmax H 1/2 2 + Ca 2 ' 1 ' 

In this expression the pumping activity of the P-Ca state which has only one calcium ion 
bound is neglected. The high cooperativity prevents us from deciding whether this state is 
contributing to the total pump activity or not. 

Equation (|THl has the form of a Hill equation with Hill coefficient n = 2. t\ can be 
identified with the maximal pump rate J m ax- Since the maximum rate of the calcium 
pumping is limited by r\ in equation dJ, r\ = J m ax stands for the saturating properties of 
the pumping process. H\i2 = is the half activation concentration and a measure for 

the affinity of the pump to calcium. We will use equation (fTljl for all further investigations 
in this paper. 

So far we have only looked at the calcium-dependent pump activity in the absence of 
calmodulin without considering pump stimulation. The activity in the unstimulated state 
corresponds to the base activity in the absence of calmodulin. It is very likely that the 
interaction with calmodulin causes the increase in pumping activity since it is well-known 
that calmodulin increases both the affinity for calcium as well as the maximal activity 
and both effects occur during stimulation Q|. The transition dynamics is given by the 
stimulation model in equations ©-(EI). Hence, the total pumping rate of a system becomes 

T ■ Cn 2 T* ■ Cn 2 

j _ f '-'max ^ » j, ^max r\ 

- /unstim • H ^2 + ^2 + ^tim " R ^ 2 + ^ , [^) 

where /unstim and / s tim are the fractions of unstimulated and stimulated pump with / uns tim = 

P+Px aTlH * . _ p *+ p x 
P ana j s tim — p 
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One input parameter of the stimulation model is the total pump concentration Pq. This 
concentration Pq is taken from Caride et al. [j| to be 0.005 fjM. Assuming a sufficiently 
large amount of pumps, the total pump concentration is irrelevant for the pumping model 
since we consider fractions of pumps in different states only. In contrast in the stimulation 
model the proportion between the total calmodulin concentration, Mo, and the total pump 
concentration, Po, is significant. If for example the amount of PMCA's exceeds the number 
of available calmodulin proteins not all pumps could be stimulated. 

For simulations the surface density of PMCA's has to be transfered into the concen- 
tration Pq. This implies the neglection of diffusion of the calcium-calmodulin complex 
to and from the pump, which is justified for rather small compartments only. For larger 
compartments a space resolving model with a local pump concentration has to be used. 

Pumping parameter determination with steady-state pump activity By setting 
all differential equations (equations (fTB|l . lfT7|) and (|T8|). see Appendix) of the stimulation 
model (reactions CO)-©) to zero the steady-state distribution of P, Px, Px an d P* is found. 
This determines the asymptotic steady-state fractions / un stim and / s ti m in equation (|To|) . 
Using this distribution we calculate the asymptotic pumping activity in dependence of 
the calcium and calmodulin concentration. Many measurements of the steady-state pump 
activity in the presence and the absence of calmodulin have been performed [H, Q, 24, 
l^ . Fitting equation (|T5|) to one of these experiments yields J max , >/ max , ^1/2 an d Hy 2 . 
In figure Q] we use the steady-state activity measurements done by Caride et al. [<]]. figure Q] 

With /unstim = 1 and /stim = in equation l|15|) we fit the steady-state activity curve 
in the absence of calmodulin which supplies J max and Hi/ 2 in equation dJ. Equation 
(|To|) is used to fit the data in the presence of calmodulin (see dashed line in figure 0J. We 
determine J max and H*, 2 with that fit, whereby the steady-state fractions /unstim and / s ti m 
are determined by the stimulation model. These fractions respect the fact that at lower 
calcium concentrations in the presence of calmodulin no stimulated pump form is available 
(/stim — 0). This can be seen by inspecting the dashed-dotted line in figure 01 which is 
plotted with the artificial assumption that only stimulated pump form (/ s tim = 1) is present 
for all calcium concentrations. The gap between this virtual line (dashed-dotted) and the 
fit of equation ifTKjl (dashed line) at low calcium concentrations is determined by the choice 
of the stimulation parameter set of table El Although we determine J max and Hy 2 with 
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that fit, changes in J max and Hy 2 only alter the slope and the saturation value whereas the 
gap at low calcium concentrations is not influenced by these changes: J max an d ^1/2 do 
not alter the low calcium concentration limit. Note that the fit of the case with calmodulin 
being present, and therefore the determination of J max and H^, 2 , is accomplished with 
the Bootstrap method using the Metropolis fit algorithm that also takes into account the 
stimulation and relaxation constant data sets. The pumping parameters obtained from this 
routine are summarized in table 0] The impact of stimulation by calmodulin on PMCA 
function, namely to increase both the affinity to calcium H\/2 and the maximum pump rate 
J m£H , can be affirmed. tabled 

The obtained parameters J max and J max refer to the investigated system, i.e. they are 
determined by the ensemble of investigated pumps. In contrast to the pump rate, the affinity 
to calcium expressed by -ffi/2 or H*^ 2 is a universal parameter since it does not depend on 
the surface density. Knowing the expression level of PMCA's within the system and the 
average protein mass one can use J max and J max to calculate the maximum unstimulated 
^single and stimulated pumping rate </ s * mg i e of a single pump respectively, which is then a 
system-independent parameter characteristic for the isoform studied. Caride et al. achieved 
with the baculovirus expression system an amount of 5% PMCA's of total membrane protein 
(private communication). With this data and with an overall average protein mass of 130 
kDa Q| we deduced a turnover rate of 10.4 Hz out of J max = 0.24/zmol/(mg • min) for the 
h2b isoform and 31.7 Hz out of J max = 0.73/imol/(mg • min) for the stimulated h4b rate. 

With the use of the pumping parameters and the stimulation rate constants (table 
12} we can calculate the calmodulin dependence of the steady-state activity (see figure 03). 
The displayed fractional steady-state activity / is a universal description, i.e. it is not figure 
dependent on the surface density. The quantitative behavior of our calculation is confirmed 
by experimental data of Penheiter et al. j^] (compare our figure [3 with figure 3 on page 
17730 in this publication). Note that we calculated this steady-state activity out of the 
available data without further fitting routines. 

Results 

Comprehensive PMCA dynamics With the knowledge of the stimulation parameters 
from table 121 and the pumping properties in tabled we are now able to simulate the time 
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dependent behavior of the calcium pumping rate of the PMCA including the stimulation. 



it Pi ti 

consecutively to different calcium concentrations. The time course of Phi produced is shown 

in figure El (a) (crosses). It is rather likely that during a single turnover of the PMCA one 

ATP is hydrolyzed as well as one calcium ion transported [l,Q, [], Q|. Therefore the rate 

of Phj production directly corresponds to the calcium pumping rate. Utilizing equation figure El 

(|15j). we can calculate and sum up the time course of Phj production. With the use of the 

universal parameter sets we adapt our simulation (full line in figure El) to this experiment 

by adjusting the unknown surface density of PMCA's. This yields an expression level of 

9.8 % PMCA's of total membrane protein, which is comparable to the quoted 5 % (see 

before). Note that we have added a constant base level rate with respect to the non-zero 

Phj production rate in the absence of calcium (Ca = O^M) in figure El (a), which has 

to be related to a different source within the investigated tissue. In fact, in the absence 

of calcium no ions can be transported, i.e. no ATP can be hydrolyzed, no Phj produced. 

Note that during the second exposure to 0.5 fjM. calcium the measured activity significantly 

decreases compared to the first time at high calcium concentration, i.e. even under the same 

experimental conditions the Ph, production rate changes during the experiment. This might 

be related to dwindling resources such as ATP or the increasing significance of ATP and 

calcium diffusion. Such effects naturally cannot be reproduced by the present model. 

After 300 seconds at 0.5 fjM. calcium, during which a PMCA fraction is stimulated, the 
concentration was decreased to 0.05 fiM calcium for different durations and raised again 
to high calcium concentration. The turnover rate at the end of the second low calcium 
exposure JqqI° is, due to the gradual decay of stimulated pump form, dependent on the 

j{second) j-(first) 

duration of the low calcium phase. The fraction of change in J0.05, °' ( ° fi 5 rst ) — mffi is shown 
in figure El (b). The temporal behavior of pump function is quantitatively reproduced. 



Discussion 

Our combined model of stimulation and pumping is able to reproduce the PMCA behavior. 
Therefore the assumptions of the stimulation model could serve as a possible explanation for 
the underlying biological steps. An essential element of the model is the separation of the 
attachment of the calmodulin-calcium complex to the pump (reaction J2J) and the stimula- 
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tion, which occurs in an additional step (reaction (J3J). The convergence of the stimulation 
constant k^^ at high calcium concentrations in figure Q provides strong evidence for that 
assumption. The translocation step in reaction © has the same characteristics. These 
steps can be seen as internal conformational changes where their velocity cannot be further 
increased by higher concentrations. In contrast to reaction (j3J which is confined by the 
total available pump concentration Pq and k^, the velocity of reaction (J2J can be arbitrarily 
increased with higher calcium or calmodulin concentrations. Similarly the saturation of the 
relaxation constant fc^stim a ^ l° w calcium concentrations calls for a rate limiting reaction 
during relaxation. Regarding the relaxation pathway in this context yields two steps which 
could meet this property, i.e. forward reaction Q and reaction ©. The Metropolis fit 
algorithm provides similar outcomes for both possibilities. In general the routine delivers 
&4 < &5, i.e. reaction (@J limits the relaxation rate. The experimental data could also be 
reproduced assuming k^ > k§, which corresponds to the irreversible relaxation being the 
limiting step. However, figure 4b in Q] provides clear evidence that the limiting step is 
related to the dissociation of calmodulin, hence step (@J) in our model. 

The separation between stimulation (reactions and (0)) and relaxation (reactions 
(J3J and (JSJ) enables us to describe the different temporal behaviors independently. This 
is required since for example the isoforms h2b and h4b display opposite stimulation and 
relaxation behavior: The h2b isoform is stimulated more quickly than h4b but relaxes more 
slowly under the same conditions. Assuming the stimulation reaction (j3J would become 
reversible, the formation and degradation of stimulated pump form P x (t) could be described 
together by the first three reactions (QJ, and Q without the use of reactions (JIJ and 
(J5J. In such a case the calcium and calmodulin dependent stimulation and relaxation 
velocity could not be adjusted autonomously which is required for the reproduction of the 
experimental data. The formal introduction of two different pathways may be interpreted 
as different underlying transformational changes. 

Comparing our parameters with those of the model proposed by Penheiter et al. in 
reveals strong agreement. These authors argue that the preferred stimulation route in 
their branched model is the binding of the X4 complex to the unstimulated (closed, in their 
terminology) conformation and the subsequent stimulation (opening) of the pump. This 
corresponds to our single stimulation step. The second route assumes a stimulated (open) 
state of the pump in the absence of the X4 complex which would be stabilized by the binding 
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of the calmodulin-calcium complex. This branch had to be introduced in Q| in order to 
describe the initial increase in fluorescence. As in our model, the conformational change 
limits the stimulation. The limiting steps of the Penheiter et al. model are determined by 
k\ and k\, where k\ describes the preferred stimulation route (here p refers to reaction rates 
in the Penheiter model). The value for k\ = 0.034 1/s is similar to our k^ (see h4b isoform 



in table ^ and Their and our mode 
the calmodulin-calcium complex \l 



confirm the known affinity of the h4b isoform for 

mm. 



Since we have based the parameter deduction on inorganic phosphate, Phj, release mea- 
surements and the Penheiter model relies on fluorescence measurements based on the bind- 
ing of calmodulin to the pump, the agreement of the parameters justifies both approaches. 
However, the reduction of the Penheiter model to the stimulation pathway ending in a 
stimulated pump with isomerized CaM protein (TA-CaM-PMCAg) might not account for 
the competing stimulation and relaxation processes of the pump. This fact together with 
the branched model ansatz could be the reason for slightly different reaction rates. 

As already mentioned, the choice of k{ (i= ±2, 3, ±4, 5) and therefore the choice of K has 
a sensitive effect on the stimulation and relaxation graphs in figure ^ and [21 Furthermore 
the steady-state graphs in figureHJare strongly influenced by these parameters. The attempt 
to reproduce all experimental data for both isoforms on the basis of our model by using a 
single value of K leads us to the following conclusion: The affinity of the PMCA for the 
X4 complex is altered during the stimulation process. The unstimulated pump P exhibits 
a lower affinity than the stimulated pump P* (compare K2 = ^-2/^2 to K4 = k^/k-^). 
For both isoforms the value of K± is thousandfold lower compared to K2 (K4 = 0.0008 
nM for h2b, K4 = 0.002 nM for h4b). Within the framework of our model the assumption 
K2 = K± fails to adequately reproduce all experimental data. In the case of the isoform h4b 
this assumption does not allow a sufficient amount of PMCA to become stimulated in the 
steady-state calculation. In principle the experimental data in figure Q] could still be fitted 
by adjusting Hy 2 and J^axi but this would contradict the statement that J max is about 
20% of J^ax f° r the h4b isoform P, IS, 2l3|. Note that this result is indirectly supported 
by the observation that the Ca 2+ /calmodulin-dependent protein kinase II (CaMKII) target 
protein has a thousandfold higher affinity to the X4 complex in the stimulated state, an 
effect denoted as calmodulin trapping }2(\ \. 

The fit results in table|H]show that apart from k<i and k^ the parameters remain relatively 
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constant for different values of K. Hence, an exact determination of K under the specific 
experimental conditions could improve the determination of k 2 and k§ but would not change 
the remaining values of the stimulation, the pumping parameters or the main conclusions. 
Note, that for all of the three assumed values, i.e. K = 0.1 , 1 and 10 l//iM 4 , the 
experimental data of both isoforms could be reproduced with equal accuracy. 

We have also investigated a different approach for PMCA relaxation. In this scenario 
only calcium dissociates from the complex in step (@J)'. The following reaction 
comprises the relaxation and the detachment of calmodulin. Though the reproduction of 
the available experimental data is ensured with this approach, there are two aspects which 
distinguish it from the present model. First, the analytical expression for the relaxation 
constant becomes independent of calmodulin due to the irreversibility of step It is, 

therefore, impossible to predict the calmodulin dependence of the relaxation rate as done 
in figure |2 Second, similar to the change of the pump affinity for the X4 complex in 
the present model, the affinity of calmodulin for calcium would be different depending on 
whether calmodulin is free (reaction £Q|) or bound to the PMCA pump (reaction (@J'). We 
can find experimental support for both scenarios. As previously quoted, CaMKII is known 
to change its affinity for the calmodulin-calcium complex j^, whereas Olwin et al. report a 
change of the dissociation constant between CaM and calcium over two orders of magnitude 
in the presence and in the absence of a target protein, i.e. in their case the rabbit skeletal 
muscle myosin light chain kinase [^J. No experimental evidence is known to us whether 
calcium or the whole X4 complex detaches first from the pump. However, since we treat 
the balance between free X4, calmodulin and calcium in quasi-steady-state our approach 
presented here is the more general one since reaction (@J makes only the statement that the 
detachment of calmodulin limits the relaxation rate but does not restrict in which order 
the components dissociate. 

In the description of the pumping behavior we have incorporated the experimental 
observations that although the translocation step involves only one calcium ion (indicated 
by the 1:1 stoichiometry of ATP to transported calcium) the pump binds two Ca 2+ ions with 
high affinity . The experimental data provide strong evidence for the assumption 

that the binding of the two calcium ions to the intracellular binding site of the pump 
is highly cooperative. The deduced pumping rate expression is in accordance with these 
observations. The high cooperativity prevents us from determining whether only P-Ca2 
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or also P-Ca can perform the transport process of one calcium ion. Our general ansatz 
could in principle account for recent findings of Guerini et al. in which the Ca 2+ -dependent 
ATPase activity is fitted by a Hill function with Hill coefficient 1 However, the fit to 
the available data (see figure 0} suggests a Hill coefficient 2. 

Our aim to characterize pump isoforms by different universal parameter sets is achieved 
for the h2b and h4b isoforms. We assume the same functionality for all isoforms, hence these 
parameters express relevant differences in the stimulation and pumping behavior. With the 
concept of providing universal parameter sets we claim that the sets are applicable to other 
simulations including hPMCA2b or hPMCA4b isoforms. 

Out of the universal parameters of the h4b isoform we can calculate an accurate predic- 
tion of the calmodulin dependence of the steady-state activity. Using the expression level 
and the overall average protein mass we calculated single pump turnover rates. Based on a 
measurement of Elwess et al. in 1997 Q| and using an expression level of 0.2 % we com- 
puted a turnover rate of 8.5 Hz for the rPMCA2b isoform. Our calculated single turnover 
rates (table Ql) are in the same range as those of Blaustein who indicated a turnover of as 

n, 

30 Hz without specifying the isoform (|4j and private communication). Adamo et al. made 
direct measurements of the turnover rate and determined it to be 33 Hz for the erythrocyte 
PMCA Q. However, the calculated turnover rates, even being in the appropriate range, 
may differ from the presented values since the expression levels of the pump can not be 
determined with much accuracy, and factors such as other proteins and lipids interacting 
with the pump in the biological membrane environment can have substantial effects. 

For comparison, the turnover rate of the Sarcoplasmic Reticulum Ca 2+ -ATPase (SERCA) 
lies within the same range. The SERCA1, SERCA2a and SERCA3 isoforms transport cal- 
cium with a turnover rate of « 10 Hz, whereas the the SERCA2b isoform exhibits a rate 
of « 5 Hz Q Note that in contrast to the PMCA the SERCA carries two calcium ions 
per pump cycle reflected in the Ca 2+ :ATP stoichiometry of 2:1 [3, 13] 

The isoform specific universal set of parameters may be applied to different scenarios 
by incorporating the system specific expression levels of PMCA. However, alterations of 
the stimulation rate constants may occur under modified experimental conditions. Since 
we consider proteins with polar binding sites we cannot exclude a change of dynamical 
properties at different ionic strengths of the experimental solutions. Corresponding experi- 
ments could elucidate the importance of such effects. Also our theoretical proposition of the 
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calmodulin dependence of the stimulation and the relaxation constant which differs from 
the linear behavior found by Caride et al. calls for further experimental investigations. 
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Appendix 

Stimulation model By using the Law of Mass Action the stimulation model ©-© is 
rewritten as a set of differential equations for concentrations 



dP k 2 P(M -P x -P* x ] 



dt l + l/{KCa A ) 

dP x k 2 P(M -P x -P* x ) 



+ k- 2 Px + h(P -P-P x -P x ), (16) 

(k- 2 + k 3 )P x , (17) 



dt 1 + l/(KCa* 

k 3 P X - k A P x + k^(P -P-P X - Py) 7" : ytr - (18) 



dP *x _ , n D r, V *,M -P X -P X 



dt ~ 6 ^ "j— . ^ , X) l + l/(KCa*) ' 

Pump conservation Pq = P + Px + P x + P* where Pq is the total pump concentration and 
calmodulin conservation Mq = M+X^+P x +P x with Mq as total calmodulin concentration 
is respected. We solve this system of coupled ordinary differential equations of first order 
with a fourth-order Runge-Kutta method. This has been implemented in a C++ program. 
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Hill fit of the stimulation constant The data from [9j can be fitted by the Hill equation 

{.max . fi A 

k^L(Ca) = )^ )4 Ca , (19) 

with the maximal stimulation constant fc™?^ and the half maximum concentration Ca^^ . 
According to (JBJ), k^ t ( m converges to ks for large calcium or calmodulin concentrations (see 



below), hence k^ confines the maximum stimulation rate and can be identified with k, 



max 
stim ' 
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This reflects the fact that equation (J3J limits the stimulation rate. By setting the right 
hand side of equation equal to ^3/2, an expression for the Ca 2+ concentration at half 
maximal stimulation is derived 



r (1/2) _ 4 / fc-2 + fc 3 /2 

Cflstim " V k 2 -K-M- (20) 
Note that the calmodulin concentration M is kept constant. 

In a similar manner, the experimental relaxation constants can be fitted by an inverse 
Hill equation. Note that we expect the dependence of the stimulation and relaxation con- 
stant on calmodulin to be fitted by a Hill equation with Hill coefficient 1. 

Limes of equation (j^J We consider the limes of k^ t ( m for X4 — ► 00 in the case of real 



and positive rate constants k 2 , k_ 2 and k, 



3- 



,(2/3) fc-2 + k 3 + k 2 X A / (k- 2 + k 3 + k 2 XA 2 
x^J^ = 2 ± V I 2 j " ^ 3 * 4 



fc_ 2 + k 3 + k 2 X A 4k 2 k 3 X 4 \ 

= x^ 2 Pf- MPw j (21) 

(2) (3) 

For fcg ti this expression diverges, whereas we find lim k^J = k 3 using the rule of de 
l'Hospital. 



References 

[1] L. Zylinska and M. Soszynski, Plasma membrane Ca 2+ -ATPase in excitable and nonex- 
citable cells, Acta Biochim. Pol. 47, 529 (2000). 

[2] J. T. Penniston and A. Enyedi, Modulation of the plasma membrane Ca 2+ pump, J. 
Membr. Biol. 165, 101 (1998). 

[3] E. Carafoli, The calcium pumping ATPase of the plasma membrane, Annu. Rev. Phys- 
iol. 53, 531 (1991). 

[4] M. Juhaszova, P. Church, M. P. Blaustein, and E. F. Stanley, Location of calcium 
transporters at presynaptic terminals, Eur. J. Neurosci. 12, 839 (2000). 

[5] D. Guerini, The significance of the isoforms of plasma membrane calcium ATPase, Cell 
Tissue Res. 292, 191 (1998). 

19 



[6] T. P. Stauffer, D. Guerini, and E. Carafoli, Tissue distribution of the four gene products 
of the plasma membrane calcium pump, J. Biol. Chem. 270, 12184 (1995). 

[7] E. Carafoli, The Ca 2+ pump of the plasma membrane, J. Biol. Chem. 267, 2115 (1992). 

[8] D. Chin and A. R. Means, Calmodulin: a prototypical calcium sensor, Trends Cell Biol. 
10, 322 (2000). 

[9] A. J. Caride, A. R. Penheiter, A. G. Filoteo, Z. Bajzer, A. Enyedi, and J. T. Penniston, 
The plasma membrane calcium pump displays memory of past calcium spikes, J. Biol. 
Chem. 276, 39797 (2001). 

[10] S. Linse, A. Helmersson, and S. Forsen, Calcium binding to calmodulin and its globular 
domains, J. Biol. Chem. 266, 8050 (1991). 

[11] W. J. Pottorf, S. P. Duckies, and J. N. Buchholz, Mechanisms of calcium buffering in 
adrenergic neurones and effects of ageing: testing the limits of homeostasis, J. Auton. 
Pharmacol. 20, 63 (2000). 

[12] A. Persechini and B. Cronk, The relationship between the free concentrations of Ca 2+ 
an Ca 2+ -calmodulin in intact cells, J. Biol. Chem. 274, 6827 (1999). 

[13] T. Meyer, P. I. Hanson, L. Stryer, and H. Schulman, Calmodulin trapping by calcium- 
calmodulin- dependent protein kinase, Science 256, 1199 (1992). 

[14] A. Persechini, K. Yano, and P. M. Stemmer, Ca 2+ binding and energy coupling in the 
calmodulin-myosin light chain kinase complex, J. Biol. Chem. 275, 4199 (2000). 

[15] N. L. Elwess, A. G. Filoteo, A. Enyedi, and J. T. Penniston, Plasma membrane Ca 2+ 
pump isoforms 2a and 2b are unusually responsive to calmodulin and Ca 2+ , J. Biol. 
Chem. 272, 17981 (1997). 

[16] A. R. Penheiter, B. Zeljko, A. G. Filoteo, R. Thorogate, K. Torok, and A. J. Caride, A 
model of the activation of plasma membrane calcium pump isoform J^b by calmodulin, 
J. Bio. Chem. 42, 12115 (2003). 

[17] A. J. Caride, A. G. Filoteo, N. L. Elwess, A. K. Verma, A. Enyedi, Z. Bajzer, and J. T. 
Penniston, The rate of activation by calmodulin of isoform 4 of the plasma membrane 



20 



Ca 2+ pump is slow ans is changed by alternative splicing, J. Biol. Chem. 274, 35227 
(1999). 

[18] A. Enyedi, A. G. Filoteo, G. Gardos, and J. T. Penniston, Calmodulin-binding do- 
mains from isozymes of the plasma membrane calcium pump have different regulatory 
properties, J. Biol. Chem. 266, 8952 (1991). 

[19] W. Ba-Thein, A. J. Caride, A. G. Filoteo, A. Enyedi, K. Paszty, C. L. Croy, and 
J. T. Penniston, Chimaeras reveal the role of the catalytic core in the activation of the 
plasma membrane Ca 2Jr pump, Biochem. J. 356, 241 (2001). 

[20] A. R. Penheiter, A. J. Caride, A. Enyedi, and J. T. Penniston, Tryptophan 1093 is 
largely responsible for the slow off rate off calmodulin from plasma membrane calcium 
pump 4b, J. Biol. Chem. 277, 17728 (2002). 

[21] A. J. Caride, A. G. Filoteo, A. R. Penheiter, K. Paszty, A. Enyedi, and J. T. Penniston, 
Delayed activation of the plasma membrane calcium pump by a sudden increase in 
Ca 2+ : fast pumps reside in fast cells, Cell Calcium 30, 49 (2001). 

[22] J. Lytton, M. Westlin, S. E. Burk, G. E. Shull, and D. H. MacLennan, Functional 
comparisions between isoforms of the sacroplasmic or endoplasmic reticulum family of 
calicum pumps, J. Biol. Chem. 267, 14483 (1992). 

[23] A. K. Verma, A. Enyedi, A. G. Filoteo, E. E. Strehler, and J. T. Penniston, Plasma 
membrane calcium pump isoform 4a has a longer calmodulin-binding domain than 4b, 
J. Biol. Chem. 271, 3713 (1996). 

[24] H. Hilfiker, D. Guerini, and E. Carafoli, Cloning and expression of the isoform 2 ot the 
human plasma membrane Ca 2+ ATPase, J. Biol. Chem. 269, 26178 (1994). 

[25] H. P. Adamo and M. E. Grimaldi, Functional consequences of relocation the c-terminal 
calmodulin-binding autoinhibitory domains of the plasma membrane Ca 2+ pump near 
the n-terminus, Biochem. J. 331, 763 (1998). 

[26] A. Hudmon and H. Schulman, Neuronal Ca 2+ /calmodulin- dependent protein kinase II: 
The role of structure and autoregulation in cellular function, Annu. Rev. Biochem. 71, 
473 (2002). 



21 



[27] B. B. Olwin, A. M. Edelman, E. G. Krebs, and D. R. Strom, Quantitation of energy 
coupling between Ca2+, calmodulin, skeletal muscle myosin light chain kinase, and 
kinase substrates, J. Biol. Chem. 259, 10949 (1984). 

[28] D. Guerini, A. Zecca-Mazza, and E. Carafoli, Single amino acid mutations in trans- 
membrane domain 5 confer to the plasme membrane Ca 2+ pump properties typical of 
the Ca 2+ pump of endo(sacro)plasmic reticulum, J. Bio. Chem. 275, 31361 (2000). 

[29] H. Adamo, A. Rega, and P. Garrahan, Pre- steady- state phosphorylation of the human 
red cell Ca2+-ATPase, J. Bio. Chem. 263, 17548 (1988). 



22 



Tables 



parameter 


h2b 


h4b 


unit 


source 


K 


> 0.029 


1/juM 4 


see text 




0.055 


0.024 


1/s 


l, max 
"'stim 


&4 


0.015 


0.0347 


1/s 


k cxp 

relax 



Table 1: Stimulation parameters and ranges determined by Hill equation fits 

Pump specific rate constants and ranges determined by fits of the Hill equation to the 
experimental data from Caride et al. (2001) 



parameter 


h2b 


c.i.(68%) 


h4b 


c.i.(68%) 


unit 




1.99 


1.92. ..8.91 


0.096 


0.088.. .0.168 


1/OuM-s) 




0.056 


0.053.. .0.058 


0.023 


0.021. ..0.025 


1/s 




0.016 


0.015. ..0.052 


0.035 


0.035. ..0.147 


1/s 


/c_4 


20001 


19988.. .20010 


19999 


19978.. .20017 


l/(/xM-s) 


k 5 


0.12 


0.02.. .0.35 


0.84 


0.14.. .0.97 


1/s 



Table 2: Determination of stimulation parameters with a Metropolis algorithm 

The medians and the confidence intervals (c.i.) of 68 % are calculated with a Bootstrap 
method combined with a Metropolis algorithm using experimental data for the stimulation 
constant, relaxation constant and steady-state pumping activity (see text for more details). 
K is chosen to be 1 /xM -4 . /c_2 can be calculated from hz using the known dissociation 
constant if 2 = fc-2/^2 (see text). 
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parameter 


h2b 


h4b 


unit 


K 


0.1 


1 


10 


0.1 


1 


10 


1///M 4 




24.8 


1.99 


0.20 


1.46 


0.096 


0.0093 


1/OM-s) 




0.058 


0.056 


0.056 


0.023 


0.023 


0.023 


1/s 




0.015 


0.016 


0.018 


0.036 


0.035 


0.035 


1/s 




19991 


20001 


20001 


19993 


19999 


19998 


1/OuM-s) 


h 


0.019 


0.12 


0.68 


0.06 


0.84 


7.4 


1/s 


J* 

"mas 


0.24 


0.24 


0.24 


0.73 


0.73 


0.76 


/imol 
mg-min 


IT* 
H l/2 


0.44 


0.45 


0.45 


0.56 


0.57 


0.54 


fjM 



Table 3: Parameters at different values of K 

See table [21 0] and the text for more details concerning the determination of these values. 



parameter 


h2b 


c.i.(%68) 


h4b 


c.i.(%68) 


unit 


characteristic 


'An ax 


0.116 (48% of J* „) 


±0.008 


0.148 (20% /j,; ax ) 


±0.018 


/imol 
mg-min 


system specific 


^single 


5.0 


±0.3 


6.4 


±0.8 


Hz 


universal 


T* 

^ max 


0.24 


0.23. ..0.26 


0.73 


0.68. ..0.75 


fimol 
mg-min 


system specific 


7* 

single 


10.4 


10.0. ..11. 3 


31.7 


29.5. ..32. 5 


Hz 


universal 


-#1/2 


0.63 


±0.07 


1.45 


±0.26 




universal 


H l/2 


0.45 


0.29. ..0.53 


0.57 


0.46. ..0.62 




universal 



Table 4: Pumping parameters 

Parameters of the steady-state pump activity determined using experimental data of Caride 
et al. 2001 Q|. K is chosen to be 1 /zM -4 . The activity values in the absence of calmodulin 
are determined by a gnuplot fit to the experimental data. The medians and the confidence 
intervals (c.i.) of J^ax an d H*, 2 are derived from a Bootstrap method combined with the 
Metropolis algorithm (see text for more details). 
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Figure legends 

(3) 

Figure 1 k st { m is plotted against calcium at constant calmodulin concentration of 0.117 
/iM for isoforms h2b and h4b according to equation (j^J (full lines). The dotted lines 

(3) 

represent the fcg tim dependence on calmodulin at fixed calcium concentration of 0.8 fiM 
for the h2b and 1 (jM. for the h4b isoform. The stimulation rate k c *? m (squares) has been 
measured at various calcium concentration with constant calmodulin concentration of 0.117 
(with kind permission of Caride [j]]). 



Figure 2 The relaxation constant & re iax is plotted for varying calcium (full lines) and 
calmodulin (dotted lines) concentrations where calmodulin (M = 0.117 fj,M for both iso- 
forms) and calcium (Ca = 0.25 /xM for h2b; Ca = 0.3 fiM for h4b) are kept constant, 
respectively. Corresponding measurements at M = 0.117 fjM. are shown (squares). With 



kind permission of Caride 



Figure 3 The dynamics of the fractions / un stim = P ~p[ x (full lines for h2b, dashed- 
dotted lines for h4b) and / s tim = — prp^ (dashed lines for h2b, dotted lines for h4b) are 
shown during stimulation (panel (a)) and relaxation (panel (b)). The total calmodulin 
concentration is 0.117 fiM and Pq = 0.005 fj,M. Note that the formation and degradation of 
stimulated and unstimulated pump form depend on the available calcium and calmodulin 
concentration. The present calcium concentration is 0.5 fjM during stimulation (a) and 0.1 
fjM during relaxation (b). 

Figure 4 Steady-state calcium dependent pump activity of both isoforms. The measure- 
ments in the absence (triangles) and in the presence (circles) of 0.117 uM. calmodulin have 

Q Q 

been kindly provided by Caride [3|. The total pump concentration Pq is taken from |2J to 
be 0.005 /iM. Without calmodulin (full line) the fit of equation (fl8|) to the data displays 
no difference to the fit of the simplified equation (fT4^1 which presumes a cooperative bind- 
ing. Therefore, only the fit of equation lfT4^) is shown. The dashed line shows the realistic 
case and displays / s tim — for low calcium concentrations. The dashed-dotted line depicts 
the pump activity with the assumption that only stimulated pump form is present at all 
calcium concentrations, i.e. / s tim = 1 and /unstim = 0. 
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Figure 5 The fractional calmodulin dependent steady-state pump activity / = (J — 
«^max)/(<^max ~~ ^max) of the h4b isoform at Ca = 0.7 fiM is shown. J max is the activity in 
the absence of calmodulin and J^ax 1S the steady-state activity in the presence of saturating 
calmodulin, both at the accordant calcium concentration of 0.7 /xM. The calculation has 
been done with P = 0.005 /jM. 

Figure 6 Figure (a) shows the theoretical (full line) and the experimentally measured 
(crosses) time course of Phj production during different calcium concentration exposures, 
shown on the top of the panel. The present calmodulin concentration was 0.117 /iM. 

(^■(second) ^j(first) \ 
- i25 — AJ 05 J as a function of time 

at low calcium (0.05 //M), experimentally measured (crosses) and theoretically predicted 

(full line). Both experiments were performed with the h2b isoform (with kind permission 

of Caride Q). 
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(b) h4b isoform 



Figure 1: 
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(b) h4b isoform 



Figure 2: 
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Figure 3: 
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(a) time course of Phi production 
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